REPORT  DOCUMENTATION  PAGE 


Form  Approved  OMB  No.  0704-0188 


Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  the  burden,  to  Department  of  Defense,  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports  (0704-0188),  1215  Jefferson 
Davis  Highway,  Suite  1204,  Arlington,  VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  penalty  for  failing  to  comply 
with  a  collection  of  information  if  it  does  not  display  a  currently  valid  OMB  control  number. 

PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 


1.  REPORT  DATE  (DD-MM-YYYY)  2.  REPORT  TYPE 

21-12-2004  Final  Report 

3.  DATES  COVERED  (From  -  To) 

4  March  2003  -  04-Jan  -05 

4.  TITLE  AND  SUBTITLE 

Detection  and  Characterization  of  Phase-Coded  Radar  Signals 

5a.  CONTRACT  NUMBER 

FA8655-03-1  -3030 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

Dr.  Ernest  R.  Adams 

5d.  PROJECT  NUMBER 

5d.  TASK  NUMBER 

5e.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Cranfield  University  (RMCS) 

Shrivenham  SN6  8LA 

United  Kingdom 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

N/A 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

EOARD 

PSC  802  BOX  14 

FPO  09499-0014 

10.  SPONSOR/MONITOR’S  ACRONYM(S) 

11.  SPONSOR/MONITOR’S  REPORT  NUMBER(S) 

SPC  03-3030 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  is  unlimited. 

13.  SUPPLEMENTARY  NOTES 

14.  ABSTRACT 

This  report  results  from  a  contract  tasking  Cranfield  University  (RMCS)  as  follows:  The  Grantee  will  investigate  various  mathematical 
concepts  of  detection  and  characterization  of  simulated  radar  signals.  As  detailed  in  the  technical  proposal,  at  least  four  approaches  will  be 
investigated:  1.  A  varied  matched  filter  and  spectral  analysis.  2.  A  neural  network  (NN)  approach,  based  on  the  use  of  NNs  to  detect  pulse 
radar.  3.  Methods  derived  from  CDMA  (code  division  multiple  access)  multiuser  detection.  4.  Independent  component  analysis  (ICA).  He  will 
employ  MATLAB  and  SIMULINK  to  generate  signals  and  MATLAB  for  their  analysis.  Efficiency  of  algorithms  will  be  an  important 
consideration  in  assessing  the  performance  of  the  mathematical  approaches  developed. 


15.  SUBJECT  TERMS 

EOARD,  Signal  Processing,  radar 


16.  SECURITY  CLASSIFICATION  OF: 


a.  REPORT 

b.  ABSTRACT 

c.  THIS  PAGE 

UNCLAS 

UNCLAS 

UNCLAS 

17.  LIMITATION  OF 
ABSTRACT 
UL 


18,  NUMBER 
OF  PAGES 

21 


19a.  NAME  OF  RESPONSIBLE  PERSON 

PAUL  LOSIEWICZ,  Ph.  D. 


19b.  TELEPHONE  NUMBER  (Include  area  code) 

+44  20  7514  4474 


Standard  Form  298  (Rev.  8/98) 

Prescribed  by  ANSI  Std.  Z39-18 


CONTRACT  TITLE: 


SPONSOR: 


SPONSOR’S  REFERENCE: 


FINAL  REPORT: 


Detection  and  Characterization  of  Phase-coded 
Radar  Signals 


AFRL/SNRP,  WPAFB  through  EOARD,  London 


033030 


Research  carried  out  up  to  December  2004  by 
Dr  E  R  Adams,  Cranfield  University  (RMCS) 


1 .  Summary 


Relatively  simple  correlation/matched  filter  techniques  were  developed  to  detect  and 
identify  phase-coded  radar  signals  based  on  true  Barker  codes  and  longer,  combined 
Barker  codes.  Detection  and  identification  rates  were  good,  though  not  so  good  as 
those  possible  with  more  sophisticated  techniques  like  neural  networks.  However,  the 
simplicity  and  computational  efficiency  of  the  new  methods  make  them  worthwhile 
for  some  applications. 

2.  Description  of  the  problem 

The  problem  is  to  detect  and  characterize  pulsed  phase-coded  radar  signals  in  noise. 
Gaussian  noise  is  assumed,  although  further  work  would  consider  non-Gaussian  noise 
and  interference  from  DS/SS-based  CDMA  transmissions.  Partial  knowledge  of  the 
radar  signal  is  assumed,  in  particular  the  carrier  frequency  and  the  type  of  pseudo¬ 
random  sequence  used  to  determine  phase.  Barker  codes  are  investigated  first  but,  as 
the  longest  known  true  Barker  code  is  only  13  bits  and  short  codes  require  high  power 
levels  for  a  given  energy,  longer  codes  created  from  true  Barker  codes  are  also 
considered.  The  extended  problem  of  detecting  multi-phase  Barker  codes,  other 
pseudo-random  codes  like  m-sequences  and  forms  of  wideband  interrupted  CW  was 
not  attempted. 
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3.  Properties  of  phase-coded  radar  signals 

A  medium/long-range  surveillance  radar  typically  requires  a  range  resolution  A R  of 
15  m.  This  AR  implies  a  waveform  RF  bandwidth  B  of  10  MHz.  From  Gabor’s 
theorem,  bit  duration  is  1  IB  =  1 0  7s.  Assuming  the  13-bit  Barker  code,  the  pulse  in 
Figure  1  results. 
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Figure  3.1:  Radar  pulse  for  a  13-bit  Barker  code 


Time  delay  td  is  related  to  range-to-target  R  by  the  equation 

2  R 

ctd=2R  =>  td= — , 


(3.1) 


so  time-delay  resolution  A  td  is  related  to  range  resolution  AR  by  the  equation: 


A  td  - 


2A R 


AR  = 


cA  t„ 


c  2 

But  A  td  is  known  to  be  MB ,  giving  range  resolution 


(3.2) 


A/i  =  Ai  =  _L 

2  2  B 


(3.3) 


Thus  5=10  MHz  gives  the  desired  range  resolution  of  15  m.  Pulse  duration  r  results 
from  the  choice  of  the  13-bit  Barker  code. 


The  pulse-repetition  interval  T,  illustrated  in  Figure  3.2,  is  a  more  fundamental 
variable  than  z.  It  is  determined  by  the  maximum  unambiguous  range,  Ru : 


T  = 


2  R„ 


c 


(3.4) 


r  should  then  be  chosen  to  satisfy  3  z<T  in  order  to  avoid  eclipsing.  In  practice,  for 
short  pulses  like  the  13-bit  Barker  code,  T»  3  r.  Greater  efficiency  results  (power 
reduced)  if  z approaches  773  (z=  T produces  CW),  which  provides  the  motivation  for 
examining  the  longer  codes  in  section  6. 
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Figure  3.2:  Train  of  regularly  repeating  pulses 

The  length  of  the  pulse  train  is  to  =  nT,  i.e.  the  train  has  n  pulses.  The  variable  n  is 
determined  by  the  requirement  that  2EINo  at  maximum  range  must  be  20  dB,  where 
E  =  Svi,  S  being  signal  power.  This  may  alternatively  be  stated  in  the  relation 
between  output  and  input  SNRs: 


2  E 


(3.5) 


i.e.  InBris  the  integration  improvement  factor. 

Finally,  Ts  is  the  interval  between  pulse  trains,  satisfying: 

r,=T-(„  (3.6) 

ilh 


where  Q  is  the  solid  angle  of  search  and  Q*  is  the  solid  angle  of  the  radar  beam.  In 
practice,  the  choice  of  Ts  is  more  arbitrary  than  for  other  parameters.  In  simulations, 
Ts  was  initially  assumed  to  be  constant.  Later,  it  was  assumed  to  be  unifonnly 
distributed  between  Ts  -  to/2  and  Ts  +  to/2. 

The  following  values  were  used  in  initial  simulations: 


B=  10  MHz, 


1 

B 


(bit  length)  = 


10^7  s 


r=1.3xl06s  (for  13 -bit  Barker  code) 


t=  m  x  10  7  s  for  pseudo-random  /?z-bit  code 


2R 

T  =  1  x  10~4  s,  based  on  Ru  =  100  km,  T  =  — - 

c 

«  =  20  =>  to  =  nT=  1.4  x  10~2  s 


Ts  =  3  s  (constant) 


Ts  =  3  +  (rand  -  0.5)to,  where  rand  is  the  Matlab  function  returning  a 
uniformly  distributed  value  between  0  and  1 . 
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4. 


Barker  and  other  pseudo-random  codes  for  phase  coding 


4.1  True  Barker  codes 

The  nine  known  Barker  codes  are  represented  in  Table  4.1. 


Length  N 

Cn 

2 

+  + 

2 

-  + 

3 

+  +  - 

4 

+  +  -  + 

4 

+  +  +  - 

5 

+  +  +  -  + 

7 

+++--+- 

11 

+  +  + - - 

13 

+++++--++-+-+ 

Table  4.1:  Barker  codes 

Barker  codes  have  autocorrelation  functions  (ACFs)  with  equal  time  sidelobes:  they 
are  the  only  codes  for  which  ACF  sidelobes  at  zero  Doppler  have  levels  <N  (code 
length).  As  explained  in  section  2,  they  are  too  short  for  some  applications  -  the 
greatest  sidelobe  reduction,  for  N=  13,  is  -22.3  dB.  However  it  is  possible  to 
combine  Barker  codes  to  create  longer  codes  with  good  properties. 

4.2  Combined  Barker  codes 

This  method  involves  phase  coding  with  one  Barker  code  within  each  segment  of 
another  Barker  code.  For  example,  the  Barker  code  of  length  3,  B3,  may  be  combined 
with  that  of  length  5,  B5,  to  create  a  code  of  length  15  (Figure  4.1): 

B53  =  (+  +  +  -+,  +  +  +  -  +, - +  -)  (4.1) 

Note  that  Bnm  is  the  Kronecker  product  Bn®Bm ,  regarding  sequences  as  row  vectors. 


Figure  4.1:  Combined  Barker  Code  B53 

The  order  in  which  the  two  codes  are  combined  is  significant,  i.e.  Bn,m  ^  Bm,n.  [1,2] 
give  examples  of  the  different  code  properties  of  the  alternative  orderings.  E.g.  B413, 
for  which  each  bit  of  the  13 -bit  code  Bi3  is  coded  into  four  bits,  has  a  zero  Doppler 
ACF  with  four  side  peaks  of  amplitude  13  at  shifts  ±1  and  ±3,  and  12  peaks  of 
amplitude  4.  Bi3(4  has  the  same  number  of  side  peaks  of  amplitude  >1,  but  the  side 
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peaks  of  amplitude  13  occur  at  shifts  of +13  and  ±39.  The  main  peak  of  the  ACF  in 
both  cases  has  amplitude  52.  B443  is  useful  if  expected  interference  is  well  separated 
in  range  from  the  target. 

Combining  pairs  of  Barker  codes  simply  in  this  way  can  give  codes  up  to 

11x13  =  143  bits  long.  The  properties  of  3-way  combinations,  e.g.  3x5x7,  are  not 

known.  Note  that  Bnn  is  not  used  because  of  high  partial  correlations. 

4.3  n-Phase  Barker  codes 

//-phase  Barker  sequences,  e.g.  using  phase  alphabets  such  as  multiples  of  n!6,  are 
defined  and  their  properties  explored  in  [3,4],  They  will  be  investigated  if  time 
allows. 
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5. 


Methods  for  detection  and  characterization  of  Barker  sequences 


5 . 1  Correlation/pattern  recognition  techniques 

These  newly  developed  techniques  begin  with  calculating  the  ACF  of  a  length  Tp  of 
signal,  where  Tp>  n  T\  the  highest  reasonable  value  of  the  unknown  nT,  and  Tp«Ts, 
i.e.  Tp  is  selected  so  there  is  a  high  probability  it  will  include  a  single  pulse  train  of  n 
pulses  of  length  T.  FFTs  are  used  to  estimate  the  ACF  to  improve  efficiency. 

Figure  5.1  illustrates  the  ACF  for  a  noiseless  pulse  train  of  B3  pulses. 


Figure  5.1:  ACF  of  3-length  code 

Pattern  recognition  algorithms  are  under  development  to  recognize  the  repeating 
waveform,  linearly  decreasing  in  amplitude  at  multiples  of  \T\,  with  the  aim  of 
estimating  T.  Alternatively,  pattern  recognition  may  be  carried  out  in  the  frequency 
domain.  Assuming  signal  x{t)  is  sampled  at  intervals  of  At,  i.e.  x,  =  x(iAt),  and 
T=  MAt,  then  the  pulse  train  of  n  Barker  codes  A,  0  <  /  <  N  -  1,  may  be  written: 


yt  =  A  *  L S(i-kM) 

k= 0 


(5.1) 


The  DFT  of  the  pulse  train  v,  is  thus  the  product  of  the  DFT  of  the  Barker  code  with 
the  DFT  of  the  finite  train  of  impulse  functions.  If  n  — >  00,  a  line  spectrum  of  the 
Barker  code  results,  with  spectral  line  spacing  \!T  Hz.  For  finite  n,  this  spectrum  is 
convolved  with  sine  ( nTf ),  i.e.  a  very  narrow  sine  function  with  zeros  at  integer 
multiples  of  1  In  T.  Assuming  T=  MAt  and  Tp  =  PAt,  the  frequency  resolution  of  the 
FFT  for  window  Tp  is  HP  At,  so  the  spectral  lines  at  1/ MAt  will  be  sampled  only  if  M 
divides  P.  This  is  illustrated  in  Figure  5.2.  Algorithms  are  being  developed  to 
estimate  T  from  the  DFT  or  power  spectrum  of  x,. 
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nMAt 


Figure  5.2:  Fourier  transform  of  signal  with  7^-length  window 

The  estimate  T  ,  either  from  the  ACF  or  FFT  method,  is  used  to  identify  the  unknown 
waveform  by  dividing  the  interval  Ta  into  K  segments  of  length  7  and  averaging  ( K  is 
the  largest  integer  for  which  Kf  <  Ta ).  Assuming  f  =  M  'At , 

ui  =  fxi+w  0<i<M-\  (5.2) 

K  k= o 

The  averaging  reduces  noise  and  allows  easier  detection  of  the  single  Barker 
waveform  in  u,  with  a  matched  filter  or  neural  network. 

5.2  Simulation  results  for  correlation  techniques 

The  performance  of  a  combination  of  the  techniques  in  5. 1  was  assessed  in  initial 
simulations.  To  reduce  computational  effort,  T  was  taken  to  be  7  x  1CT5  s.  Otherwise, 
the  values  in  section  3  were  used.  A  two-stage  method  detected  the  presence  of  a 
repeating  pulse-waveform  and  estimated  T,  the  pulse-repetition  interval.  Because  of 
the  high  cost  of  ACF  estimation,  detection  and  initial  estimation  of  T  were  based  on 
the  signal’s  Fourier  spectrum.  This  allowed  a  more  precise  estimate  of  T,  following 
detection  and  initial  T  estimation,  from  a  narrow  search  of  ACF  values.  The  refined  T 
estimate  was  then  used  to  produce  an  averaged  pulse-wavefonn,  i.e.  one  with  greatly 
reduced  noise.  The  final  stage  was  to  identify  the  waveform  by  correlation  with 
Barker  codes  -  a  matched  filter  approach.  In  more  detail,  the  stages  were: 
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(1)  As  described  in  5. 1,  the  DFT  of  a  signal  containing  a  repeating  waveform 
consists  of  narrow  sine  functions  spaced  at  l/T  Hz  intervals,  the  sine  amplitudes 
depending  on  the  waveform’s  Fourier  transform.  The  FFT  At  of  the  F’-lcngth  signal 
(same  symbols  used  as  in  5.2)  was  used  to  test  for  repeating  sine  functions  by  finding 
the  value  i  =  r  which  maximizes 


So 


Y;  = 


Xi  +  X2i  +  X3i  +  '  • '  +  XIV 

k ■ 


(5.3) 


i*r  =>  Yr>Yi  (5.4) 

This  Yr  was  then  tested  against  a  threshold  to  detect  the  presence  of  a  repeating 
waveform.  If  a  repeating  waveform  was  detected  the  following  approximate  estimate 
of  T  was  calculated: 


PAt 

r 


(5.5) 


The  estimate  is  not  exact  as,  in  general,  1  IT  is  not  a  sampled  frequency  of  the  FFT,  i.e. 
the  sine  functions’  centres  are  not  sampled. 

Detection  rates  for  Bn  pulse-trains  in  alternative  noise  environments  are  given  in 
table  5.2.  False  alann  rates,  i.e.  detection  rates  when  no  pulse-train  is  present,  are  also 
given. 

(2)  Given  T{  from  the  first  stage  (5.5),  a  more  accurate  estimate  of  T  (and  one  that 
may  be  exact)  may  be  derived  by  examining  ACF  values  for  shifts  close  to  Tx .  For 
such  a  corresponding  m,  a  shift  index  for  which  mAt  =  Jj ,  the  ACF  is 

Cm=  E[AA+m],  (5-6) 


estimated  by  averaging  XiXi+m  over  all  possible  i  ( 1  <i  <P  -  m).  For  m  values 
satisfying 


mAt  -  Tl 


<  AT, 


(5.7) 


AT  being  detennined  by  trial  and  error,  the  following  statistic  was  maximized: 


Dm  = 


Cm  +C2m  +Ci,m  +---  +  C, 


kmm 


(5.8) 


For  M'  corresponding  to  this  maximum,  i.e. 


m*M'  ^  Dm,  >  Dm 


(5.9) 


the  final  estimate  for  T  becomes: 


T  =  M'At  (5.10) 

(3)  T  was  then  used  to  produce  a  waveform  with  reduced  noise  by  averaging 
signal  segments  of  length  T  (5.2).  The  waveform  was  identified  by  correlating  the 
averaged  segment  with  alternative  Barker  codes.  A  match  was  found  if  the  maximum 
correlation  exceeded  a  pre-detennined  threshold.  Although  only  B 13  signals  were 
simulated,  the  averaged  segment  was  also  correlated  with  Bn  and  false  detection  rates 
(deciding  Bn  was  present  instead  of  the  actual  B 13)  were  estimated.  Also  estimated 
were  failure  rates  to  detect  any  pulse-train  when  B 13  was  present.  The  identification 
results  may  be  seen  in  table  5.3.  The  confusion  between  B 13  and  Bn  may  be 
explained  by  their  high  cross-correlation: 


CCj  (1 3,1 1)  =  £  i?13  (i)Bl  j  (i  +  j) 

i 

=  0---0 -10-1-2-1  0-3-410-3730-1-4301-2101  0---0 

(5.11) 

More  sophisticated  pattern  recognition  algorithms  are  being  developed  to  distinguish 
Barker  codes.  The  cross-correlation  properties  of  Barker  codes  are  also  under 
investigation. 


Noise 

%  Detection 

%  False  alarms 

6  dB 

96 

2 

3  dB 

86 

8 

0  dB 

62 

42 

Table  5.2:  Detection/false  alarm  rates  for  Bn 


Noise 

%  Identification  (Bn) 

%  False  identification  (Bn) 

%  Failure  to 
identify 

6  dB 

100 

4 

0 

3  dB 

98 

12 

0 

0  dB 

90 

26 

6 

Table  5.3:  Identification  rates  for  Bn 


Each  estimated  %  in  tables  5.2  and  5.3  is  based  on  50  simulations.  Note  that  a  0  dB 
noise  level  corresponds  to  noise  with  the  same  mean  amplitude  (i.e.  1)  as  the  Barker 
code.  However,  the  noise  persists  throughout  the  signal  at  this  level.  3  dB 
corresponds  to  a  mean  noise  amplitude  of  0.707  (power  0.5)  and  6  dB  to  0.5  (power 
0.25) 
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Note  also  that  a  radar  receiver  needs  an  output  SNR  (i.e.  2E/No)  of  100  (20  dB),  and 


So  the  receiver  needs 


5  _  100 
N  n(Br) 


(5.12) 


(5.13) 


In  the  above  case,  n  =  20  and  Bz=  13,  so  the  required  S/N=  0.38,  i.e.  a  mean- 
amplitude  ratio  of  0.62. 
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6.  Methods  for  detection  and  identification  of  combined  Barker  sequences 
6.1  Correlation/matched-filter  techniques 

Combined  Barker  sequences  may  be  expressed  as  the  Kronecker  product  of  a  pair  of 
Barker  sequences.  This  is  worthwhile  as  the  Matlab  routine  ‘kron’  may  be  used  to 
produce  combined  sequences.  The  Barker  sequences  of  lengths  r  and  5  may  be 
represented  by  the  vectors: 

Br  =(«,■)=  («t  «2-«r) 

b,  =(/*,)=(fl &-A)  ' 

A  combined  Barker  sequence  B,s  is  produced  by  the  Kronecker  product  Br  0  Bv ,  i.e. 
taking  Bs  =  (fij)  and  replacing  each  entry  fy  by  the  r- vector  /?/«/): 


B,v  =  B,®B,  =  (rk) 


=  (AK)AK)-^SK)) 

(6.2) 

=  (P\<*\  P\ai  ■■■P\Ocr  P2<*\  P1a1...p2ar.  .  . 

Psa  1  Psa2  ■  ■  ■  Psar  ) 

Correlating  the  combined  Barker  sequence  Brv  with  the 
the  matched-filter  output: 

inner  sequence  Br  produces 

1  v1 

ck  ~  j+\Y k+j 

r  j= o 

(6.3) 

Note  that  Ck  produces  fi,  values  every  r-th  value  of  k: 

Cy  —  — Py{(Xy  +  CC  +  .  .  .  +  CC  ^ 
r 

cl+r  =~P2[a{  +  a2  +  ...  +  «;) 
r 

(6.4) 

ClHs-l)r  =~fis{ai+a 2  +  •  •  •  +  ar  ) 
V 

As,  Vi,  at  =  ±1  =>  a,2  =1,  all  the  (. ..)  terms  in  (6. 4)  equal  r,  giving: 

C\+ jr  =  P j+\  5  0  <  y  <  j  - 1  (6.5) 

For  all  other  Ck  terms: 
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Thus,  correlating  every  r-th  value  of  (c*)  with  the  outer  sequence  BiS  produces  the 
maximum  matched-filter  output: 


1  s-i 

cc{  =  -  X/?/+,cl+/, 

S  j= 0 


=  —  X/?/  ,  from  6.5 

■s  j=i 

=  1  ,  as  p)  =  1 


(6.7) 


The  above  suggests  a  2-dimensional  search  for  the  presence  of  B„  in  a  sample. 
Successively  searching  for  B,.  and  Bs  is  more  efficient  than  directly  searching  for  B,v, 
i.e.  for  the  presence  of  a  repeating  (%)  in  the  sample.  Correlating  with  an  r-length 
sequence  and  an  5-length  is  clearly  more  efficient  than  correlating  with  one  of  length 
rs.  Furthermore,  it  is  possible  to  estimate  r  before  going  on  to  estimate  s.  Because  of 
poor  cross-correlation  properties  of  some  Barker  sequences,  initial  correlation  of  the 
sample  with  an  incorrect  Br  may  lead  to  a  weak  repeating  pattern,  increasing  the 
likelihood  of  false  identification. 

6.1.1.  Estimation  of  B j  components 

For  a  particular  B,.  and  sample  T=(vi),  the  following  matched  filter  results: 


1  ?'-l 

ck  =~Tocj+lyk+j  ,  k>  1  (6.8) 

r  j= o 


The  first  stage  in  estimating  the  outer  sequence  is  to  correlate  alternative  Bs  with  c/c 
values  at  k-intervals  of  r ,  producing  the  cck  sequence  (matched-filter  output): 


cel  =  -  I Pj+\Ck+jr  ,  k  >  1  (6.9) 

s  j= o 

If  the  correct  BiS  is  used,  the  maximum  value  of  1  is  produced  when  that  Bv  is  aligned 
with  identical  r-spaced  c  values,  as  in  (6.7),  which  happens  once  for  each  Bra 
waveform  (every  M-th  value  of  k  in  (6.9),  T  =  MAt).  Thus,  for  noiseless  (y,): 


max[cc/)  ]  =  1 

k 


and,  for  other  k' : 


1  ck'+jr  1^  -  and  Pj  =  ±1  => 

cck, 

= 

1  S_1 

P j+\C k'+ jr 

r 

s  j= 0 

(6.10) 


(6.11) 
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Tests  for  Bs  exploit  the  expected  pulse  train  of  1  ’s  in  ccsk  for  the  correct  s.  The  train 
of  n  pulses  at  intervals  of  T  (=MAt),  beginning  at  mT,  may  be  represented: 


m+n-l  .  x  it  —  tn  1  00  / 

I  S(t  -  /T)  =  rect  l  —  >  Xf>(t 

i—m  I  I  /=— on 


iT) 


Its  Fourier  transform  is: 


-y'2^o  <■ 


r 


=  zze  J  v'”smc(nTf  )*  X  dj  / - 

i=— oo  \  T 


Its  amplitude  spectrum  is  thus: 


r 


|FT{  J  =  n|sinc(«r/)|*  'LS\f~- 

i=- oo  y  1 


(6.12) 


(6.13) 


(6.14) 


i.e.  a  series  of  absolute  sine  functions  centred  on  frequencies  0  Hz,  ±  —  Hz,  ±  —  Hz 

1  2 

etc.  The  sine  centred  on  0  Hz  has  zero  values  at  frequencies  ±  — ,  ±  —  Hz  etc, 

nT  nT 

where  —  «  —  .  Tests  for  Bs  may  thus  be  based  on  the  amplitude  spectrum  of  cc). 
nT  T 

-  a  phase  independent  test,  i.e.  ^-independent. 


An  alternative  and  simpler  test  estimates  M  by  taking  average  of  ccsk  values  at 
^-intervals  of  k.  Different  values  for  p  are  tried  within  the  likely  range  for  M,  which 
may  be  derived  by  examining  the  amplitude  spectrum  of  ccsk  as  above,  or  may  be 
known  a  priori.  The  /^-interval  averages  are: 

1  lp~1  N 

ccs-  =  —  Z  cc-  -  ,  L  <  (6.15) 

Lp  y=o  ^  p 


where  N  is  the  number  of  available  y  samples.  The  largest  of  the  p  possible  values 
(z  =  1,2 ,...,p)  is  taken  for  each  possible p  and  the  overall  maximum  over p  is 
compared  with  a  threshold  to  detect  5: 

max  [ccsf  ]  >  threshold  ?  (6.16) 

p,\<i<p 
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A  final  alternative  is  to  test  the  series  resulting  from  taking  every  r-th  value  of  the  ck 
in  (6.8)  using  the  techniques  of  section  5,  i.e.  test  for  a  true  Barker  code. 

Tests  for  a  single  true  Barker  code  or  a  combined  Barker  code  in  the  sample  may  be 
combined.  Using  an  initial  B,-  matched  filter  may  result  in  a  pulse  train  of  1  ’s  or  a 
recurring  pattern  associated  with  a  further  Bs.  Thus  it  is  possible  to  decide  there  is 
just  a  repeating  true  Barker  code  Br  present  or  go  on  to  detect  the  additional  Bs  of  a 
combined  code. 

6. 1 .2.  Effects  of  additive  Gaussian  noise 

Applying  the  matched  filter  of  (6.8)  to  a  sample  consisting  of  a  signal  (x,)  in  Gaussian 
noise  («,): 


yi=xi+ni,  (6.17) 

where  are  independent  samples  from  a  zero-mean  Gaussian  distribution  of  variance 
cr,  produces 


1  ''-1  /  \ 
ck  =  -  I aj+\\?k+j+nk+j) 

r  j= o 

1  r- 1  1  r—1 

^jOtj+l^k+j  ^  j+lM/c+j 

r  j= o  r  7=o 


(6.18) 


2 

As  aj+ 1  =  ±1,  the  noise  in  ck  has  zero-mean  and  variance  a  !r.  This  effect  is  repeated 
for  cc‘l  in  (6.9),  which  has  a  variance  of  o2/rs.  This  clear  aids  detection,  particularly 

for  longer  sequences,  as  does  the  smaller  value  for  cc  when  /( is  not  perfectly  aligned 

(6.11). 

6.2  Simulation  results  for  matched  filter/correlation  techniques 


As  an  illustration,  identification  rates  are  given  for  B53  and  Bnj.  ck  and  ccsk  were 

calculated  as  in  (6.8,  6.9)  and  tests  for  Bl5  used  (6.15).  M  was  assumed  to  be  in  the 
range  600-800,  when  the  actual  value  was  700. 


Noise 

%  Identification 

%  False  identification 

%  Failure  to 
identify 

6dB 

84 

13 

11 

3  dB 

81 

14 

16 

OdB 

70 

23 

27 

Table  6.1:  Identification  rates  for  B5,3 
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Noise 

%  Identification 

%  False  identification 

%  Failure  to 
identify 

6dB 

91 

17 

7 

3  dB 

88 

21 

10 

OdB 

79 

34 

21 

Table  6.2:  Identification  rates  for  Bi3,7 

Each  estimated  %  in  tables  6.1  and  6.2  is  based  on  100  simulations,  with  the  same 
noise  assumptions  as  described  in  5.2. 

Clearly,  identification  is  better  for  the  longer  B^j,  as  expected,  but  false 
identification  is  also  higher  due  to  higher  cross-correlations  for  the  constituent  Barker 
codes. 
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7. 


Alternative  methods  for  Barker  sequence  detection 


7. 1  Neural  network  detection 

Neural  networks  have  previously  been  trained  to  recognize  particular  codes  [5].  The 
current  investigation  has  focussed  on  a  3-layer  perceptron  network  with  13  input 
elements  trained  to  detect  B13.  The  output  layer  has  one  neuron,  which  may  be  shown 
to  give  a  sufficient  statistic.  All  neurons  in  the  hidden  and  output  layers  use  a  sigmoid 
function,  and  neuron  thresholds  are  zero. 


Sequence 

Training  input 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

2 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

3 

-1 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

4 

1 

-1 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

5 

-1 

1 

-1 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

13 

1 

1 

1 

1 

-1 

-1 

1 

1 

-1 

1 

-1 

1 

0 

14 

1 

1 

1 

1 

1 

-1 

-1 

1 

1 

-1 

1 

-1 

1 

15 

0 

1 

1 

1 

1 

1 

-1 

-1 

1 

1 

-1 

1 

-1 

24 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

1 

1 

25 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

1 

26 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

Table  7.1:  Training  set,  without  noise,  for  Bi3 

Training  was  carried  out  with  the  26  possible  sequences  in  Table  7. 1  with  additive 
noise,  to  attempt  to  get  output  1  for  sequence  14  and  0  for  all  other  sequences. 

[5]  has  shown  that  after  training  with  noiseless  codes,  the  signal-to-sidelobe  ratio  of 
this  kind  of  network  can  achieve  42.7  dB  for  Bn  and  49.7  dB  for  the  63-bit 
m-sequence,  much  higher  than  the  respective  17  dB  and  22  dB  for  matched  filters. 
The  neural  network  also  performed  much  better  in  the  presence  of  non-Gaussian 
noise.  Current  results  so  far  for  B13  with  additive  Gaussian  noise  have  achieved 
around  35  dB,  still  a  much  superior  detector  to  the  conventional  matched  filter. 

Future  work  will  investigate  the  effects  of  the  number  of  neurons  in  the  hidden  layer 
and  alternative  learning  algorithms.  Clearly,  different  networks  would  be  needed  for 
different  codes,  and  subsequent  analysis  required  to  detennine  other  parameters, 
e.g.  T(or  n).  An  obvious  disadvantage  of  neural  networks  is  their  computational 
complexity. 
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7.2  Methods  derived  from  CDMA  multiuser  detection 


The  pulse-compression  technique  in  phase-coded  radar  involves  the  transmission  of  a 
long-duration  wide-bandwidth  signal  code,  properties  shared  by  CDMS  transmissions. 
Thus  more  severe  problems  arise  when  interfering  CDMA  signals  are  present.  In 
these  cases,  it  may  be  beneficial  to  use  a  modified  multiuser  detector  (MUD), 
designed  to  reject  multiple  access  interference  (MAI)  in  CDMA,  possibly  blindly  and 
adaptively.  MUDs  use  information  about  all  users’  codes  to  cancel  MAI,  i.e.  all  users’ 
signals  except  the  wanted  user’s:  in  general,  the  more  information  used  the  better,  and 
slower,  the  detection.  In  this  application,  all  users  constitute  the  MAI  and  the  radar 
signal  the  desired  user.  MUDs  using  subspace  methods  [6]  will  be  investigated:  the 
basis  of  this  approach  is  that,  by  perfonning  an  eigendecomposition  of  the 
autocorrelation  matrix  of  signal  samples,  the  signal  is  separated  into  radar  signal  and 
CDMA  user  subspace  and  the  noise  subspace. 

MUDs  exploiting  the  cyclostationarity  and  resulting  high  degree  of  spectral 
correlation  of  signals  will  also  be  investigated.  Silent  periods  between  pulse  trains 
may  be  exploited  to  derive  the  (cyclostationary)  CDMA  interference  model,  removing 
the  need  for  training  sequences  [7]. 

7.3  Independent  component  analysis  (ICA) 

This  is  a  statistical  technique  which  is  suitable  for  blind  estimation  and  detection 
when  signals  from  multiple  sensors  are  available.  ICA  finds  a  linear  transfonnation  of 
the  signal  vector  to  separate  it  into  its  component,  statistically  independent,  parts. 
Phase-coded  radar  signals,  any  CDMA  users’  signals  present,  and  noise  will  be 
mutually  independent.  The  following  model  is  assumed  for  the  signal  vector  y: 

y  =  Mx  +  n  (7.1) 

where  x  is  a  vector  containing  the  radar  signal  and  CDMA  users’  signals,  M  is  a 
mixing  matrix,  and  n  an  additive  noise  vector.  ICA  computes  a  linear  transform  W 
such  that  the  elements  of 


5  =  Wy  (7.2) 

are  as  statistically  independent  as  possible:  these  5-elements  are  estimates  of  the 
independent  source  signals,  i.e.  the  radar  and  CDMA  users’  signals. 

ICA  has  been  used  as  a  MUD  for  CDMA  with  little  knowledge  of  even  the  wanted 
signal  [8],  and  also  for  separating  monopulse  radar  signals  [9]. 
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8. 


Conclusions  and  future  work 


Only  the  simple,  low-complexity  correlation  and  matched-filter  methods  of  sections  5 
and  6  have  been  investigated  in  detail.  They  produce  good  detection  and 
identification  rates  for  true  Barker  codes,  less  good  for  combined  Barker  codes. 
Further  investigation  into  the  cross-correlation  properties  of  true  Barker  codes  is 
needed  to  derive  theoretical  perfonnance  limits.  It  seems  unlikely  that  the  slightly 
improved  performance  of  more  sophisticated  methods  like  neural  networks  would 
justify  the  extra  complexity.  Neural  networks  are  not  real-time  methods  for  this 
application. 

The  techniques  of  section  6  may  in  principle  be  extended  to  analyse  signals 
employing  Barker-like  codes  constructed  from  three  true  Barker  codes,  B„®BV®BW. 
E.g  Bi3®(Bn®B7)  and  (Bi3®Bn)®B7  are  codes  of  length  1001.  Fiowever,  the 
properties  of  such  codes  are  unknown  and  need  investigation. 

Finally,  alternative  noise  environments  need  to  be  explored,  in  particular  modified 
MAI  techniques  to  cancel  interfering  signals. 
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10.  Glossary 


Abbreviations 

ACF 

autocorrelation  function 

CDMA 

code  division  multiple  access 

CW 

continuous  wave 

DFT 

discrete  Fourier  transform 

FFT 

fast  Fourier  transform 

ICA 

independent  component  analysis 

MAI 

multiple-access  interference 

m-sequence 

maximal-length  sequence 

MUD 

multi-user  detector 

RF 

radio  frequency 

Symbols 

B 

Bn 

Bn,m 

C 

Cm 

CCj 

E 

m 

n 

N 

No 

R 

A R 

Ru 

s 

S 

to 

td 

A  td 

T 

Ts 

T 

Q 

Qb 


0 


bandwidth 

Barker  code  of  length  n 

combined  Barker  code  from  B„  and  Bm 

velocity  of  light 

ACF  value  for  shift  m 

Cross-correlation  function  value  for  shift  j 

input  energy,  i.e.  St0  (S  is  the  mean  input  signal  power) 

metres 

number  of  pulses  in  train 

noise  input  to  receiver,  i.e.  NqB 

input-noise  power  spectral  density  (one-sided) 

range  to  target 

range  resolution 

maximum  unambiguous  range 

seconds 

input  signal  power 
length  of  pulse  train 
time  delay 

time-delay  resolution 
pulse-repetition  interval 
interval  between  pulse  trains 
pulse  duration 
solid  angle  of  search 
solid  angle  of  radar  beam 
approximately  equal  to 
convolution 

Kronecker  product,  e.g.  (ai  a2)  0  (bi  b2)  =  (biai  bia2  b2ai  b2a2) 
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